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Multi-hadron operators are crucial for reliably extracting the masses of excited states lying above 
multi-hadron thresholds in lattice QCD Monte Carlo calculations. The construction of multi-hadron 
operators with significant coupling to the lowest-lying states of interest involves combining single 
hadron operators of various momenta. The design and implementation of large sets of spatially- 
extended single-hadron operators of definite momentum and their combinations into two-hadron 
operators are described. The single hadron operators are all assemblages of gauge-covariantly- 
displaced, smeared quark fields. Group-theoretical projections onto the irreducible representations 
of the symmetry group of a cubic spatial lattice are used in all isospin channels. Tests of these 
operators on 24^ x 128 and 32"^ x 256 anisotropic lattices using a stochastic method of treating 
the low-lying modes of quark propagation which exploits Laplacian Heaviside quark-field smearing 
are presented. The method provides reliable estimates of all needed correlations, even those that 
are particularly difficult to compute, such as rjrj — > rjr) in the scalar channel, which involves the 
subtraction of a large vacuum expectation value. A new glueball operator is introduced, and the 
evaluation of the mixing of this glueball operator with a quark-antiquark operator, tttt, and rjr] 
operators is shown to be feasible. 
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I. INTRODUCTION 



Markov-chain Monte Carlo estimates of quantum chro- 
modynamics (QCD) path integrals defined on a space- 
time lattice are a promising means of calculating the mass 
spectrum of excited-state hadron resonances. Because of 
the way in which stationary-state energies are extracted 
from the temporal correlations of suitable quantum field 
operators in such calculations, the energy of a partic- 
ular state of interest can only be determined after the 
energies of all lower-lying states in the same symmetry 
channel are obtained. Multi-hadron states populate the 
spectrum below most of the resonances of interest. To 
reliably determine the energies of such states, the use 
of appropriate multi-hadron operators is crucial. Multi- 
hadron operators with significant coupling to the low- 
lying states of interest can be obtained by combining 
single-hadron operators of various momenta. The con- 
struction and testing of single-hadron operators of defi- 
nite momentum and their combinations into two-hadron 
operators are the subject of this work. 

Our approach to constructing single baryon operators 
of zero momentum was previously described in Ref. l|. 
A slightly different method was reported in Ref. 3- 
Our first study of the nucleon and A excitations in the 
quenched approximation was presented in Ref. Q, and 
nucleon results for two flavors of dynamical quarks ap- 
peared in Ref. Q. A survey of excited-state energies in 
small volume for the isovector mesons and kaons using 



Nf — 2 + 1 dynamical quarks was given in Ref. [s] , along 
with results for the A, S, S baryons. To extend our efforts 
into larger volumes and towards u, d quark masses yield- 
ing lighter pions, the issue of multi-hadron states was ad- 
dressed in Ref. p . A new stochastic method of treating 
the low-lying modes of quark propagation which exploits 
Laplacian Heaviside quark-field smearing was presented 
in that work, although the method was briefly introduced 
with preliminary testing in Refs. [1, 0, Q- Other recent 
progress in calculating excited-state energies in lattice 
QCD can be found in Rcfs. [iH. 

In this work, the approach of Ref. fi| is extended to me- 
son operators of zero momentum and to both meson and 
baryon operators having definite nonzero momentum. A 
new glueball operator is also introduced and tested. To 
simplify our spectrum calculations as much as possible 
and to increase the statistical precision of our results, we 
make use of single-hadron operators that transform irre- 
ducibly under all symmetries of a three-dimensional cubic 
lattice with periodic boundary conditions. Our method 
of constructing such operators is described in detail in 
this paper. Spectrum results obtained using these op- 
erators will be presented in later reports, although we 
present some testing of these operators here. Our op- 
erator design utilizes group-theoretical projections. The 
point and space groups we use are well known, and the 
properties of their irreducible representations are widely 
available in the literature. However, we collect together 
and present in this paper some of the specific group the- 
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ory details needed for our operator construction for the 
convenience of the reader and as a record of our conven- 
tions and notation. 

The Monte Carlo method commonly employed in QCD 
computations applies only to space-time lattices of finite 
extent. Hence, our goal is to obtain the stationary-state 
energies of QCD in a cubic box using periodic boundary 
conditions. In such a cubic box, we no longer have full 
rotational symmetry, even in the continuous space-time 
limit. The stationary states cannot be labelled by the 
usual spin- J quantum numbers. Instead, the stationary 
states in a box with periodic boundary conditions must 
be labelled by the irreducible representations (irreps) of 
the cubic space group, even in the continuum limit. 

This paper is organized as follows. Our approach 
to building single-hadron operators of definite momen- 
tum is described in Sec. |lTl We construct operators 
that transform irreducibly under all symmetries of a 
three-dimensional cubic lattice. Monte Carlo calcula- 
tions are required to test our operator construction, and 
we present the implementation details of such compu- 
tations in Sec. IIIII Tests of our single-hadron operators 
using the stochastic LapH method on 24^ x 128 and 
32'^ X 256 anisotropic lattices with pion masses 171^^ ss 390 
and 240 MeV are then presented in Sec. IIVI How we 
combine these operators to form two-hadron operators is 
described in Sec. |Vl and initial tests are presented. In 
Sec. IVIl a new glueball operator is introduced and we 
demonstrate the feasibility of evaluating the mixing of 
this glueball operator with a quark-antiquark operator, 
an rjr] operator, and multiple two-pion operators. Con- 
cluding remarks are given in Sec. I VIII along with our 
plans for future work. 

II. SINGLE-HADRON OPERATORS OF 
DEFINITE MOMENTUM 

We extract the finite-volume stationary-state ener- 
gies of QCD from matrices of temporal correlations 
Cijitp - to) = {0\TOiitF)Oj{to) |0), where T denotes 
time-ordering, the source operators Oj(io) create the 
states of interest at an initial time tg, and the sink opera- 
tors Oi{tp) annihilate the states of interest at a later time 
tp. The correlation functions C'ij{t) can be expressed in 
terms of "path" integrals over quark fields ip, ip and gluon 
fields U involving the QCD action having the form 

S[^,^P,U]^^K[U]^P + Sg[U], (1) 

where K[U] is known as the Dirac matrix and >S'g[[/] is 
the gauge-field action. We use an anisotropic space-time 
lattice in which the temporal spacing at is smaller than 
the spacing as in the three spatial dimensions. Each 
time slice is a three-dimensional cubic lattice. We as- 
sume periodic boundary conditions in all three spatial 
directions. In order to estimate the path integrals using 
the Monte Carlo method, it is necessary to work in the 
imaginary time formalism. This also has the advantage 



of converting oscillatory factors in real time into decay- 
ing exponentials. Our Euclidean space-time conventions, 
including the Euclidean Dirac 7-matrices, are given in 
Ref. In the four- vector a;^ describing a given lat- 
tice site, the /i = 4 component specifies the position in 
time, and the fi = 1,2,3 components specify the position 
along the Cartesian x, y, z spatial directions. As usual in 
lattice gauge theory, the gluon field is introduced using 
the parallel transporter U^{x) given by the path-ordered 
exponential of the gauge field along a link in the /z di- 
rection connecting neighboring sites of the lattice. The 
Dirac spinor field "0^^ {x) annihilates a quark and creates 
an antiquark at lattice site x, where A refers to the quark 

flavor, a refers to color, and a is the Dirac spin index, 
^ 

and the field i.^) annihilates an antiquark and creates 
a quark. Unlike in Minkowski space-time, and ip must 
be treated as independent fields. When used in path in- 
tegrals, the link variables are S'C/(3) matrices and ip, 
are complex Grassmann fields. 

Our hadron operators are constructed using spatially- 
smoothed link variables Uj{x) and spatially-smeared 
quark fields ^{x). The spatial links are smeared using 
the stout-link procedure described in Ref. [IJ] . Note that 
only spatial staples are used in the link smoothening; no 
temporal staples are used. The smeared quark field for 
each quark flavor is defined by 

'tfaaix) ^ Sab{x,y) iJbaiy), (2) 

where x, y are lattice sites, a, b are color indices, and a 
is a Dirac spin component. We use the Laplacian Heav- 
iside (LapH) quark-field smearing scheme introduced in 
Ref. [23] and defined by 

5 = 9 (a^ + a) , (3) 

where A is the three-dimensional gauge-covariant Lapla- 
cian defined in terms of the stout-smeared gauge field 
U, and (Ts is the smearing cutoff parameter. More details 
concerning this smearing scheme are described in Ref. [Q] . 

All of our single-hadron operators are assemblages 
of basic building blocks which are gauge-covariantly- 
displaced, LapH-smeared quark fields: 

qLj = D^'^i^ii^ , qtj = ^aa 74 , (4) 

where a is a color index, a is a Dirac spin component, A 
is a quark fiavor, 74 is the temporal Dirac 7-matrix, and 
D^-^^ is a gauge-covariant displacement of type j. The 
displacement type is a sequence of p spatial directions on 
the lattice j — {ji,j2, ■ ■ ■ ,jp)- This displacement can be 
trivial (j = meaning no displacement), a displacement 
in a given single spatial direction on the lattice by some 
number of links (typically two or three), or a combination 
of two or more spatial lattice directions. If we define 
dr = ii + J2 + • ■ ■ + jr-i: then the displacement Z?'^' is 
defined as a product of smeared link variables: 

D^^'>(x,x') = Ujjx) Uj^{x+d2) Uj,{x+d3)... 
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FIG. 1: The axes Ox, Oy, Oz, Oa, Ob, Oc, Od, Oe, Of, 
Oa, op, O7, 05 corresponding to the group elements Cnj of 
the octahedral group O, where Cnj denotes a proper rotation 
through angle 2n /n about axis Oj. 



(5) 



The use of 74 in Eq. ((4]) is convenient for obtaining baryon 
correlation matrices that are Hermitian. To simplify no- 
tation, the Dirac spin component and the displacement 
type are sometimes combined into a single index in what 
follows. 

We can simplify our spectrum calculations as much 
as possible by working with single-hadron operators that 
transform irreducibly under all symmetries of a three- 
dimensional cubic lattice of infinite extent or finite ex- 
tent with periodic boundary conditions. Such symme- 
tries form the simple cubic space group known as 0\ in 
Schonflies notation or Pm3m in international notation. 
This crystallographic space group is a semi-direct prod- 
uct of the abelian group of allowed translations on a sim- 
ple cubic lattice and the orthogonal point group Oh- For 
bosonic systems with zero strangeness, we add G-parity 
as a symmetry operation. 

In order to construct such hadron operators, we first 
need to identify all of the symmetry operations and to 
determine how our basic building blocks transform un- 
der these operations. An element of the space group 0\ 
is here denoted by {R, b), where i? is a spatial rotation or 
reflection about the origin and is an element of the point 
group Oh, and b is an allowed spatial shift. The group 
element [R, b) corresponds to the coordinate transforma- 
tion X Rx + b. The covariantly-displaced quark fields 
transform according to 

UiRM)qtA^)U(RM) = S{R)-lq%R,{Rx + b), (6) 
t^(i?.,6)9aaj(a;)C/(Vb) - laPB^j{R^ + b)S{R)po., (7) 

where f/(_R,b) denotes the quantum field operator that 
effects the transformation [R, b), and the transformation 
matrices for spatial inversion Is and proper rotations Cnj 
through angle 27r/n about axis Oj are given by 



S{Cnj) = exp 

S{Is) = 74, 



(8) 
(9) 



with ujki = —'^TTSjki/n and uj^k = LOki = (w^jy is an 
antisymmetric tensor which parametrizes rotations and 
boosts). A rotation by tt/2 about the y-axis is conven- 
tionally denoted by C^y, and C^z denotes a rotation by 
7r/2 about the z-axis. These particular group elements 
are given by 



S{C,y) = -L(l+^,^3), 



S{C,z) = -i=(l- 



-7271)- (10) 



The allowed rotations on a three-dimensional spatially- 
isotropic cubic lattice form the octahedral group O which 
has 24 elements. For the convenience of the reader, the 
rotation axes corresponding to these group elements are 
shown in Fig. [TJ Inclusion of spatial inversion Is yields 
the point group Oh which has 48 elements occurring in 
ten conjugacy classes. All elements of Oh can be gen- 
erated from appropriate products of only C^y, 6*4^, and 
Is- Under G-parity, our basic building blocks transform 
according to 



Ug gi,(x) 4 



q%,{x){^2)po.G^^, (11) 
{12U CA^) G^^, (12) 



using the Dirac-Pauli representation for the 7-matrices, 
and where the only nonzero elements of the G flavor ma- 
trix are G"'' = -G''" = -G'*'* = 1. 

The construction of irreducible representations (irreps) 
of begins with the irreps of the abelian subgroup of 
lattice translations. These are characterized by a definite 
three-momentum p as allowed by the periodic boundary 
conditions. 

Each of our meson operators which destroys a three- 
momentum p is a linear superposition of gauge-invariant 
quark-antiquark elemental operators of the form 

as 

(13) 

where q, q are defined in Eq. (|4]) , da, dp are the spa- 
tial displacements of the q, q fields, respectively, from 
X, A, B indicate flavor, and a, (3 are compound indices 
incorporating both spin and quark-displacement types. 
The phase factor involving the quark-antiquark displace- 
ments is needed to ensure proper transformation prop- 
erties under G-parity for arbitrary displacement types. 
The "barred" operators which create a momentum p then 
take the form 



it) 



E 



e^P<-+i('i^+'i,))5,,q^^{x,t)q^a{^,t). 



Each meson sink operator has the form 

„(0 .f.AB( 



(14) 



(15) 



(or is a flavor combination of the above form), where I 
is a compound index comprised of a three-momentum p, 
an irreducible representation A of the little group of p 
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TABLE I: The spatial arrangements of the quark-antiquark 
meson operators (left) and the three-quark baryon operators 
(right). In the illustrations, the smeared quarks fields are de- 
picted by solid circles, each hollow circle indicates a smeared 
"barred" antiquark field, the solid line segments indicate co- 
variant displacements, and each hollow box indicates the lo- 
cation of a Levi-Civita color coupling. For simplicity, all dis- 
placements have the same length in an operator. 



Meson configurations 



Baryon configurations 



single-site 

O— • 

singly-displaced 



y'-displc 

u 

-displac 



doubly-displaced-L 



triply-displaced-U 



triply-displaced- O 



single-site 

singly-displaced 

•-11— • 

doubly-displaced-I 

T 

I!-* 

doubly-displaced-L 

triply-displaced-T 

triply-displaced-O 



(discussed below), the row A of the irrep, total isospin /, 
isospin projection I3, strangeness S, and an identifier la- 
beling the different operators in each symmetry channel. 
Here, we focus on mesons containing only u, d, s quarks. 
The corresponding source operators are 



AB , 



(16) 



Each of our baryon operators destroying a three- 
momentum p is a linear superposition of gauge-invariant 
elemental three-quark operators of the form 



£abc qta (x, t) Qbi3 (x, t) q^^{x, t) . 

(17) 

The "barred" three-quark elemental operators which cre- 
ate a momentum p have the form 



ABC 



Sabc q%{x,t) q"p{x,t) qaa{x,t) 



TzA 



Our baryon sink operators, being linear superpositions of 
the three-quark elemental operators, have the form 



(19) 



where again, the I label includes the momentum p, the lit- 
tle group irrep A, the row A of the irrep, isospin /, isospin 
projection /s, strangeness S, and an identifier specifying 



the different operators in each symmetry channel. We 
focus on baryons containing only u, d, s quarks. The cor- 
responding source operators are 



Bi{t) 



(20) 



In order to build up the necessary orbital and ra- 
dial structures expected in the hadron excitations, we 
use a variety of spatially-extended configurations for our 
hadron operators, as shown in Table ID First, consider the 
zero-momentum operators. The simplest meson opera- 
tors combine the quark and antiquark fields at the same 
lattice site. We refer to these as single-site (SS) opera- 
tors. In the singly-displaced (SD) meson operators, the 
quark is displaced from the antiquark along a direction 
parallel to one of the axes of the lattice. If the quark is co- 
variantly displaced from the antiquark along an L-shaped 
path, we refer to this as a doubly-displaced-L (DDL) 
operator. Displacement of the quark along a U-shaped 
path or in three orthogonal directions from the antiquark 
leads to triply-displaced-U (TDU) and triply-displaced- 
O (TDO) meson operators, respectively, as shown in Ta- 
ble m The simplest baryon operators combine the three 
quarks at a single lattice site. In the singly-displaced 
baryons, one of the quarks is displaced from the other 
two along a direction parallel to one of the axes of the 
lattice. Displacement of two quarks away from the site 
of the Levi-Civita coupling leads to doubly-displaced-I 
(DDI) and doubly-displaced-L baryon operators of zero 
momentum. All three quarks can be displaced from the 
color-coupling site, producing triply-displaced-T (TDT) 
and triply-displaced-O baryon operators of zero momen- 
tum, as illustrated in Table ID For simplicity, all displace- 
ment lengths along each of the different directions are 
taken to be the same in any given operator. For mesons, 
we use a length of Sag, and for baryons, the length is 2as- 

For nonzero momenta, we restrict our attention to 
on-axis momenta, such as in the ±x, ±y, ±£ direc- 
tions, to momenta in a planar-diagonal direction, such 
as ±x ± y, zLx ± z, ±y ± z, and momenta in a cubic- 
diagonal direction, such as ±x±y±z. We expect that the 
above momentum directions are sufficient for studying 
the stationary-states in the range of energies of interest 
to us. For on-axis momenta, we construct single-site me- 
son and baryon operators, longitudinally-singly-displaced 
(LSD) operators in which one quark is displaced along 
the direction of the momentum, and transverse-singly- 
displaced (TSD) operators, in which one quark is dis- 
placed along a direction of the lattice transverse to the 
momentum. For planar-diagonal momenta, we use single- 
site operators, transverse-singly-displaced operators in 
which one quark is displaced along the direction of the 
lattice that is perpendicular to the plane containing the 
momentum direction, and planar-singly-displaced (PSD) 
operators in which one quark is displaced along one of the 
two directions of the lattice coinciding with the nonzero 
components of the momentum. For cubic-diagonal mo- 
menta, we use single-site and singly-displaced (SD) con- 
figurations. For such momenta, displacements along the 
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TABLE II: Our choices for the reference momentum p^cf di- 
rections and the reference rotations i?ref for each momentum 
p direction that we use. 



Pref direction 


p direction 


Kc 


(0,0,1) 


( 0, 0,-1) 


C2x 




( 1, 0, 0) 






(-1, 0, 0) 






( 0, -1, 0) 


Cax 




( 0, 1, 0) 




(0, 1, 1) 


(0,-1,-1) 






(0, 1,-1) 






(0, -1, 1) 






( 1,0, 1) 


42 




(-1,0,-1) 


C2b = C2xC4z 




( 1,0,-1) 


C*2a ~ C2yC4z 




(-1,0, 1) 


Caz 




( 1, 1,0) 


C^y 




(-1,-1,0) 


C*2d = C2zC4y 




( 1,-1,0) 


C2c = C4yC2z 




(-1, 1,0) 


^4y 


( 1, 1, 1) 


( 1, 1,-1) 


Ciy 




( 1,-1, 1) 


Cax 




( 1,-1,-1) 


C2x 




(-1, 1, 1) 


Ciz 




(-1,1,-1) 


C2y 




(-1,-1, 1) 


C2z 




(-1,-1,-1) 


C2d = C2zC4y 



lattice axes are neither entirely transverse nor entirely 
longitudinal to the momentum. 

For a given flavor structure, the next step in our 
single-hadron operator construction is to find coefficients 
in Eqs. (jTS)) and (|19l) that produce operators which 
transform irrcducibly under all symmetries of the three- 
dimensional cubic lattice. First, for each class of mo- 
menta, such as on-axis or planar-diagonal, we choose 
one representative reference momentum direction Prcf- 
We then find coefficients corresponding to operators that 
transform irreducibly under the little group of Prof • Re- 
call that the little group of Pref is the subset of symmetry 
operations that leave the reference momentum p^cf h> 
variant. Next, for each momentum direction p in a class 
of momenta, we select one reference rotation Rf^^ that 
transforms Piof into p. As long as the selected group el- 
ement transforms Prcf into p, it does not matter which 
group element is chosen, but a choice must be made and 
clearly specified. Hadron operators having a momentum 
in the direction of p are then obtained by applying the 
reference rotation to the operators constructed using the 
momentum in the direction of Picf • Our choices of ref- 
erence momenta directions and reference rotations are 
listed in Table HH 

Our choices of reference directions and rotations, as 



TABLE III: The little group corresponding to reference mo- 
mentum direction (0,0,0) is Oh- The elements of the double 
group are listed below, grouped into conjugacy classes. E 
is the identity element, E represents a rotation by 2n about 
any axis, and G = EG for each element G in Oh- Spatial 
inversion 7s in Oj^ satisfies Is — E and I^^ — Is. Conjugacy 
classes Cg through Cie are not listed below. The elements of 
class Cn+8 are obtained by multiplying each of the elements 
in class d by Is- 

Ci = {E} 

C2 — { Cia , C3/3 , Csj , C3S , C'g J , C^p , Cg^ , C^l } 
C3 = {C2x,C2y,C2z,C2x,C2y,C2z} 
Ci ~ {Cix, Ciy jCiz, C^^ , , '-^4z } 

C5 = {C2a, C2b, C2C, C2d, C2e, (^2/, 

C2a, C26, C2c, C2d, C2e, ^2/} 

Ce^{E} _ _ _ _ _ _ 

C7 = {Csa, Csfj, C:jj, C35, , Ca^ , C^ry , C^s } 

Cs = {Cix, Ciy, Ciz, , Ciy , Ciz } 



well as our choices of the irreducible representation matri- 
ces, described later, are dictated mainly by simplicity. An 
alternate approach would be to choose the z-direction as 
the single reference direction, obtain all other momenta 
using a rotation defined by the Jacob- Wick convention, 
and use irreducible representation matrices correspond- 
ing to helicity states. Since our procedure for combining 
the single-hadron operators into multi-hadron operators 
is automated using Maple, we found that there was no 
great advantage in using the standard Jacob- Wick con- 
vention with helicities. 

The little groups associated with our choices of refer- 
ence momentum directions are listed in Tables HID IIV( 
IVl and IVII To describe both mesons and baryons, we 
need the single-valued and double-valued (spinorial) ir- 
reps of these groups. The double- valued representations 
of a group G arc constructed by extending the group el- 
ements to form the so-called "double group" ■ This 
is done by introducing a new generator, denoted by E, 
which represents a rotation by 27r about any axis. For 
each element G of the original group, the double group 
contains another element G — EG. For the convenience 
of the reader, the elements of the double groups associ- 
ated with our choices of reference momentum directions 
are explicitly listed in Tables Hill IIVI IVl and lVIi grouped 
into their conjugacy classes. 

The irreducible representations of these little groups 
and their characters are listed in Tables IVIIl I Villi IIXI 
and|Xl One-dimensional single- valued irreps are labelled 
by A or B, two-dimensional irreps are denoted by E, 
and three-dimensional irreps are labelled by T. One- 
dimensional double- valued irreps are denoted by F, two- 
dimensional spinor irreps are denoted by G, and four- 
dimensional irreps are indicated by -ff. A subscript g 
indicates an even-parity irrep, whereas a subscript u in- 
dicates an odd-parity irrep. For mesons which are eigen- 
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TABLE IV: The little group corresponding to reference mo- 
mentum direction (0, 0, 1) is Ca^. The elements of the double 
group C2u for this reference momentum direction are listed 
below, grouped into conjugacy classes. E is the identity ele- 
ment, E represents a rotation by 2-k about any axis. 

Ci^{E} _ 

C'Z ~ {C2z, C22} 
Cs ~ {C4z, C^^} 

C4 ~ {Is C2x , Is C2y , IsC2x, IsC2y} 
C5 — {IsC2a, IsC2b, IsC2a, IsC2b} 

C&^{E} 

C7 = {Ciz, Ciz } 



TABLE V: The little group corresponding to reference mo- 
mentum direction (0, 1, 1) is Czv The elements of the double 
group C2V for this reference momentum direction are listed 
below, grouped into conjugacy classes. 

Ci = {£} _ 

C2 = {C2e, C2e} 

C3 — {IsC2f,IsC2f} 

Ci — {IsC2x, IsC2x} 

Cs = {E} 



states of C?-parity, a superscript indicates an irrep 
describing states even under G-parity, and a "— " super- 
script indicates an irrep associated with states odd under 
G-parity. Our notation differs from that of Refs. [26l.[27|. 

Our method of constructing the hadron operators that 
transform irreducibly under each little group makes use 
of group-theoretical projections and is described in detail 
in Ref. The first step in the method is to identify a 
basis of hadron elemental operators {t) at a single time 
t that transform into one another under the elements of 
the httle group Q. The key formula in obtaining the linear 
combinations Oi of these basis operators that transform 
irreducibly under Q is 

Ot\t) ^^Y. (^) Ur Mt) Ul (21) 
flee" 

where C/^ is the double group of G, R denotes an element 

TABLE VL The little group corresponding to reference mo- 
mentum direction (1, 1, 1) is C^v The elements of the double 
group ClJ, are listed below, grouped into conjugacy classes. 

Ci = {E} 

C2 = {C35,C^/} 

C3 = {IsC2b, IsC2d, IsCzf} 

Ci = {E} 

Cs = {Csa, } 

Ce = {IsC2b, IsCzd, IsCzf} 



TABLE VIL Characters x oi the single-valued and double- 
valued irreducible representations A of the group Oh- Only 
the even-parity irreps (subscript g) and classes Ci to Cg are 
shown below. For the even-parity irreps, Xn+s — Xny where 
Xn denotes the character of A for all group elements in class 
C„. For the odd-parity irreps (subscript u instead oi g), Xn" = 





"or n — 


1 ... 8, and Xn" 


= — 


Xn' for 


n = 9. 


..16. 




A 


xt 




x^ 


X^ 


X^ 


X^ 


X^ 


x^ 




1 


1 


1 


1 


1 


1 


1 


1 


A23 


1 


1 


1 


~1 


-1 


1 


1 


-1 




2 


-1 


2 








2 


-1 





Tuj 


3 





-1 


1 


-1 


3 





1 


T2, 


3 





-1 


-1 


1 


3 





-1 


Gig 


2 


1 





^/2 





-2 


-1 




G2g 


2 


1 










-2 


-1 


V2 


Hg 


4 


-1 











-4 


1 






TABLE VIII: Characters x'^ for the single- valued and double- 
valued irreps A of the group Civ 



A 


X^ 


X^ 


X^ 


xt 


xt 


x^ 


X^ 


Ai 


1 


1 


1 


1 


1 


1 


1 


A2 


1 


1 


1 


-1 


-1 


1 


1 


Bi 


1 


1 


-1 


1 


-1 


1 


-1 


B2 


1 


1 


-1 


-1 


1 


1 


-1 


E 


2 


-2 











2 





Gi 


2 





V2 








-2 


-V2 


G2 


2 





-V2 








-2 


V2 



of Q , ggn is the number of elements in Q , (Ia is the 
dimension of the A irreducible representation, and r^(i?) 
is the matrix for clement R in irrep A. 

To carry out the projections in Eq. (|2ip. explicit repre- 
sentation matrices F are needed for each group element, 
and a representation for the Dirac 7-matrices must be 
selected. We use the Dirac-Pauli representation for the 
7-matrices, as described in Ref. [1]. Often, helicity states 
are used for moving hadrons in continuous space-time. 
We could find no great advantage in using a helicity rep- 
resentation for our subsequent lattice calculations involv- 
ing multi-hadron operators. The use of helicity states 
does not result in any reduction in the number of ir- 
rep rows that must be evaluated for the single moving 
hadron states. All of our group theory manipulations are 
implemented using Maple, so the computational effort 
is independent of our choices of representation matrices. 
Instead, we chose the simplest possible matrices for our 
irreducible representations. Our choices of representa- 
tion matrices are summarized in Tables IXIl IXIIl IXIIIi 
and IXIVI The representation matrices for all group ele- 
ments can be obtained by suitable multiplications of the 
matrices shown in these tables. 

For a given set of elemental hadron operators that 
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TABLE IX: Characters "xf^ for the single-valued and double- 
valued irreps A of the group C2v ■ 



A 




X^ 


X3 


x^ 




^1 


1 


1 


1 


1 


1 


^2 


1 


1 


-1 


-1 


1 


Si 


1 


-1 


1 


-1 


1 




1 


-1 


-1 


1 


1 


G 


2 











-2 



TABLE X: Characters x for the single-valued and double- 
valued irreps A of the group C3,, . 



A 


xt 


X^ 


X^ 


X^ 


X^ 


X^ 


Ai 


1 


1 


1 


1 


1 


1 


A2 


1 


1 


-1 


1 


1 


-1 


E 


2 


-1 





2 


-1 





Fi 


1 


-1 


i 


-1 


1 


—i 


F2 


1 


-1 


—i 


-1 


1 


i 


G 


2 


1 





-2 


-1 






TABLE XI: Our choices for the representation matrices F of 
the single-valued and double- valued irreps A of the group Oh 
for zero momentum operators. Only the even-parity irreps 
(subscript g) are shown below. The matrices for C4y and Caz 
for the odd-parity irreps (subscript u) are the same. Spatial 
inversion I,, is the third generator. r^^\ls) is the identity 
matrix for the even-parity irreps, and minus one times the 
identity for the odd-parity irreps. The matrices for all other 
group elements can be obtained from appropriate multiplica- 
tions of the matrices below and the matrix for 7.. 



A 



r(^)(C4 



A2g 



tig 



'2g 



[1] 
[-1] 

1 Vs 
-1 

001" 
010 
-1 0. 

0-1 
0-1 
100 



1 

V2 



1 -1 
1 1 



1 

75 



[1] 
[-1] 
-1 
1 

) -1 

) 1. 

1 

1 
0-1 

1-i 
1+i 



T2g 



-1 
71 



-1 
71 



Ha 



1 



1 -Vs Vs 
Vs -1 -1 

V3 1-1 

. 1 Vs Vs 



-1 
1 



1 

7^ 



1-i 
1+i 

1-i 
1-i 

l+i 

-1+i 



TABLE XII: Our choices for the representation matrices F of 
the single- valued and double- valued irreps A of the little group 
C41, for momentum in the direction (0, 0, 1). The matrices for 
all other group elements can be obtained from appropriate 
multiplications of the matrices below. 



A 



F'^) (C4 



F''^) {IsC2y 



Ai 
A2 

Bi 
B2 



[1] 

[1] 
[-1] 
[-1] 

-1 

1 



1 


1-i 


71 


. 1 + i 


-1 


' 1-i 


71 


1 + i 



[1] 
[-1] 
[1] 
[-1] 

1 
-1 

-1 

1 

-1 

1 



TABLE XIII: Our choices for the representation matrices F of 
the single- valued and double- valued irreps A of the little group 
C21, for momentum in the direction (0, 1, 1). The matrices for 
all other group elements can be obtained from appropriate 
multiplications of the matrices below. C2e = C2zCax is a 
rotation about (0, 1, 1), and C2/ = C2yCAx is a rotation about 
(0,1,-1). 



A 


F(^) (C2e) 


r^^^ {isC2f) 


^1 




[1] 




[1] 


A2 




[1] 




[-1] 


Bi 




[-1] 




[1] 


B2 




[-1] 




[-1] 


G 


1 


-i -1 


1 


i -1 


71 


1 i 


71 


1 -i 



TABLE XIV: Our choices for the representation matrices F of 
the single- valued and double- valued irreps A of the little group 
Csv for momentum in the direction (1, 1, 1). The matrices for 
all other group elements can be obtained from appropriate 
multiplications of the matrices below. Css = CAyCiz is a 
rotation about (1, 1, 1), and C26 = C2xCaz is a rotation about 
(1,-1,0). 



A 

A2 
E 

Fi 
F2 



F<^> (C3, 



F<^> iLC2b) 



[1] 
[1] 

-1 

-Vs -1 

[-1] 
[-1] 

1-i -1- 
1-i 1+i 



1 

7i 



[1] 
[-1] 
-1 
1 

ii\ 
l-i] 

1-i 
-1-i 
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transform among one another, many of the projections 
in Eq. (I2ip vanish or lead to linearly-dependent opera- 
tors, so the final step in the operator construction is to 
choose suitable linear combinations of the projected op- 
erators to obtain a final set of independent single-hadron 
operators. These linear combinations are obtained using 
a Gram-Schmidt procedure as described in Ref. [l|. 

At the end of this entire procedure, we obtain single- 
hadron annihilation operators Bij^^^ characterized by 
total isospin /, the projection of the total isospin /a, 
strangeness 5, momentum p, little group irrep A, and ir- 
rep row A. Here, we use the index i to indicate all other 
quantum numbers and identifying information. All of 
these single-hadron operators constructed as described 
above transform under a group element (i?, b) of the 
space group 0\, in which x — Rx + b, according to 

(22) 

where the Wigner rotation is given by 

Rl^ = {R^Sr' R (23) 

and is an element of the little group of p. Note that 
the above equations apply even when R refers to spatial 
inversion Is- Eqs. (|22|) and (|23|) play a crucial role when 
forming the multi-hadron operators. The behaviors of 
our operators under G-parity and isospin rotations are 
considered below. 

In Ref. [H, the odd-parity zero- momentum baryons 
were constructed from their even-parity partner op- 
erators utilizing a particular transformation involv- 
ing charge-conjugation. For a given even-parity zero- 
momentum baryon operator Bf{t), an odd-parity opera- 
tor Bfit) can be defined by rotating the three Dirac in- 
dices using the 72 matrix and replacing the expansion co- 
efficients by their complex conjugates. This yields partic- 
ular relationships between the temporal correlation ma- 
trices of the even-parity baryons and the time-reversed 
correlation matrices of the odd-parity baryon operators, 
allowing averaging over forward and backward temporal 
propagations in some cases for increased statistics. How- 
ever, for the baryon operators having nonzero momen- 
tum, parity is no longer a good quantum number since 
it reverses the three-momentum. For this reason, we do 
not bother to apply a generalization of the above proce- 
dure in constructing the odd-parity baryon operators of 
nonzero momentum, especially since the increased sta- 
tistical precision can be obtained in other more efficient 
ways when using the stochastic LapH method. 

In constructing our light pion and kaon operators, we 
take symmetry under time reversal into account. Al- 
though the lattices we use are rather large in temporal 
extent, temporal wrap-around effects can still come into 
play for the light pious and kaons in certain situations 
where high precision is needed, such as in studying the tttt 



scattering phase shifts. Energies can be extracted with 
increased statistical precision if meson operators whose 
temporal correlations are symmetric under time reversal 
are used. Symmetry under time reversal helps in ex- 
tracting meson energies whenever temporal wrap-around 
effects become non-negligible since the functional forms 
used for fitting the data have fewer parameters, leading 
to more precise energy estimates. For meson operators, 
the energies of the states propagating backwards in time 
are the same as those of the states traveling forwards in 
time, but the couplings of a given meson operator to the 
forward-propagating states can differ from the couplings 
to the backward-propagating states if symmetry under 
time reversal is not taken into account. 

We can reduce the number of needed fit parameters if 
we use meson operators whose temporal correlation ma- 
trices satisfy C^^ (t) = C^^ {Nt — t) , where Nt is the tem- 
poral extent of the lattice, assuming periodic boundary 
conditions in time. This can be achieved if each me- 
son operator itself satisfies Miit) = 77 Mi{Nt — t), with 
= 1. Under a certain symmetry of the lattice ac- 
tion that involves time reversal, the covariantly-displaced 
LapH-smeared quark fields transform according to 

laajix) > (7475)a^i q^PjiTx), 

qtji^) qa^jiTx) {J4l5)^a, (24) 

where {Tx)j — Xj and {Tx)^ — Nt — x^. Since wrap- 
around effects are a potential problem only for the light- 
est mesons, we decided to modify only the symmetry 
channels containing the lightest pseudoscalars {Tr,K,ri). 
In each case, we projected the operators into the even and 
odd operators under the above transformation. Numeri- 
cal tests using a small number of configurations showed 
that the odd operators performed somewhat better, but 
the difference was not very significant. Hence, we dis- 
carded the even operators, and kept the odd operators. 
Note that the operator V'TsV' is odd under the above 
transformation. 

As in Ref. we work in the approximation that the 
masses of the u and d quarks are equal, and we neglect 
electromagnetic interactions. In this approximation, the 
theory has an exact isotopic spin symmetry, and states 
are characterized by total isospin /, its projection I3 onto 
a given axis, and strangeness S. Again, we consider only 
the u, d, s quarks here. Incorporating this isospin sym- 
metry into our operators is straightforward and has been 
described in Ref. [l| . The flavor structure of the hadron 
operators we use are summarized in Table IXVI Due to 
an approximate SU{3) uc?s-flavor symmetry, quark fla- 
vor combinations in meson and baryon operators are of- 
ten chosen according to the irreducible representations 
of SU{3) flavor. Such combinations are simply linear 
superpositions of the operators presented in Table IXVI 
Since we plan to obtain Monte Carlo estimates of the 
complete correlation matrices of operators including all 
allowed flavor combinations, the use of linear superposi- 
tions which transform irreducibly under SU{3) flavor is 
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TABLE XV: Flavor structure of the elemental hadron anni- 
hiliation operators we use. Each is characterized by isospin 
/, maximal I3 = /, strangeness 5", and G-parity, where appli- 
cable, expressed in terms of the gauge-invariant three-quark 
and quark-antiquark operators defined in Eqs. (|17|l and (|13|l , 
respectively. Ug is the quantum operator that effects a G- 
parity transformation. Recall that f,f',h,h',b,a are even- 
parity mesons, while ri,ri' ,u}, (j>, p,n are odd-parity mesons. 
Also, a flavored meson is named K* if its total spin J and 
parity P are both odd or both even, otherwise it is named K. 



Hadron 


I = h 


s 


G 


Annihilation operators 


A++ 


3 
2 







^uuu 


E+ 


1 


-1 




^uus 


N+ 


1 

2 







f^uud _ ^duu 

^ ap'y ^ a.p'^f 


=■0 


1 
2 


-2 




^ssu 


A° 





-1 




ap'y ap'y 


n- 





-3 






f,f',v,v' 








1 


^Tp+^ii+UGm''p+<i%)u^G 

'i-^^.+UGK^U^G 


h, h' ,uj,(j} 








-1 


^Tp+^ii-UGi'frp+^iiM 


b+,p+ 


1 





1 


+ UgKWg 


a+,7v+ 


1 





-1 




K+,K*+ 


1 

2 


1 




^ ap 




1 

2 


-1 




^ ap 



unnecessary. For example, we construct isoscalar meson 
operators having flavor content uu + dd separately from 
those having flavor content ss. Our correlation matri- 
ces then allow mixings of these operators. Our choices 
of operators described in Table IXVI are dictated by com- 
putational simplicity. In summary, we construct single- 
hadron operators that transform under an isospin rota- 
tion Rr according to 

Ur^b'XMuI = B%{t) D'-P^.^iR,), (25) 

where D^^\Rr) are the familiar Wiener rotation matri- 
ces. Our meson operators M are constructed such that 
they transform under G-parity according to 

UGMl'li{t)ul, = VKM'Jl-'it), 

UGM'Jlf,it)U^a = VAM^AXi^it), (26) 

where ?7a = 1 if >S = ±1, and when the strangeness S* = 0, 
then r]\ — ±1 depending on the G-parity superscript of 
the irrep A. 

To associate our finite-box energies with observed 
hadrons, it is necessary to know which spin- J irreps of 
the continuous group of rotations occur in which irreps 
of the octahedral point group. For the convenience of 
the reader, the spin contents of the Oh single- valued and 



TABLE XVL Continuum limit spin identification; the num- 
ber n;( of times that the A single- valued irrep of the octahedral 
point group Oh occurs in the (reducible) subduction of the in- 
teger J irrep of 5(7(2). The numbers for Aiu, A2U, Eu, Tiu, 
T2u are the same as for A\g, ^29, Eg, Tig, T2g, respectively. 



J 






-fig 


Tu 


T2C, 


U 


1 
i 


U 


U 


U 


U 


1 











1 





2 








1 





1 


3 





1 





1 


1 


4 


1 





1 


1 


1 


5 








1 


2 


1 


6 


1 


1 


1 


1 


2 


7 





1 


1 


2 


2 


8 


1 





2 


2 


2 


9 


1 


1 


1 


3 


2 



TABLE XVIL Continuum limit spin identification: the num- 
ber of times that the A double-valued irrep of the oc- 
tahedral point group Oh occurs in the (reducible) subduc- 
tion of the half-integral J irrep of SU{2). The numbers for 



Giu, G2U, Hu 


are the 


same 


as for Gig 


, G2g , Hg , 


respectively. 




^G2q 




J 






1 
1 

i 1 




1 
1 



1 
1 
1 


9 
2 

11 
2 

13 
2 

15 
2 


1 
1 
1 
1 


2 

1 2 

2 2 
1 3 



double- valued irreps are listed in Tables IXVII and IXVIIl 
respectively. These tables list the number of times that 
each irrep A of Oh appears in various J irreps of SU{2) 
subduced to the double group of Oh- Table IXVIIIl is use- 
ful for identifying which hadrons of nonzero momentum 
appear in which irreps of the little groups. The decom- 
positions of the subduced representations of Oh into the 
irreps of the little groups C^y, C^y, and C2V are given 
in this table. Table IXIXI lists the irreps of Oh in which 
various common hadrons at rest appear. 



III. IMPLEMENTATION DETAILS 

In order to test the effectiveness of the single-hadron 
operators that we have designed, Monte Carlo calcula- 
tions must be carried out. Details on how the tempo- 
ral correlations of hadron operators are evaluated us- 
ing Monte Carlo integration with the stochastic LapH 
method are presented in this section. Many of these de- 
tails have already been described in Ref. Q, so only de- 
tails in addition to those in Ref. [1] are presented here, 
along with some reiteration of important run parameters. 

Our computations use the so-called stochastic LapH 
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TABLE XVIII: Subduction ^ of the irreducible representa- 
tions of Oh to the irreducible representations of the little 
groups Civ, Csv, and C2v. 



A{Oh) 


iCiv 










Al 


Al 




Al 


Axu 


A2 


A2 




A2 


Ai„ 


Bi 


A2 




B2 


Alu 


B2 


Al 




Bi 


En 


Al © Bi 


E 




Al © B2 


Eu 


A2 © B2 


E 




A2 © Bi 


Ti, 


A2®E 


A2®E 


A2 


© Bi © B2 


Tiu 


Ai^E 


Ai^E 


Al 


© Bi © 


T2g 


B2®E 


Ai®E 


Al 


© A2 © Bi 




Bi(BE 


A2(S E 


Al 


© A2 © B2 


Glg/u 


Gi 


G 




G 


G2g/u 


G2 


G 




G 


Hg/u 


Gi ©G2 


Fi © Fa © G 




2G 



TABLE XIX: The irreducible representations of Oh in which 
various commonly-known hadrons at rest occur. 



Hadron Irrep 



Hadron Irrep 



Hadron Irrep 



P 

ao 
61 



Ai^ 
Til 
A + 

r+ 
-'19 



G 



19 



K 

LO,(j) 

fo 

Ki 

A,H 



Aiu 

Ti^ 
A + 

Tig 

Gig 



K* 
hi 



At 
Tiu 

T.-g 



method 6] and are done in a sequence of steps: (a) gener- 
ation of gauge-field configurations using the Monte Carlo 
method; (b) computation of quark sinks for various noises 
and dilution projectors using the configurations from the 
first step; (c) computation of the single meson and sin- 
gle baryon sources and sinks using the quark sources and 
sinks from the second step; (d) evaluation of the cor- 
relators using the single-hadron sources and sinks and 
their combinations into multi-hadron sources and sinks; 
and (e) analysis of the correlators to extract the energies. 
These steps are discussed below. 

A description of step (a) is given in Ref. [1^. We are 
currently focusing on three Monte Carlo ensembles: (A) a 
set of 412 gauge-field configurations on a large 32'^ x 256 
anisotropic lattice with a pion mass ~ 240 MeV, 

(B) an ensemble of 551 configurations on an 24"^ x 128 
anisotropic lattice with a pion mass TOtt 390 MeV, and 

(C) an ensemble of 584 configurations on an 24"^ x 128 
anisotropic lattice with a pion mass m^^ ^ 240 MeV. 
We refer to these ensembles as the (32^1 240), (24^1 390), 
and (24'^|240) ensembles, respectively. These ensembles 
were generated using the Rational Hybrid Monte Carlo 
(RHMC) algorithm [231, which is a Metropolis method 
with a sophisticated means of proposing a global change 
to the gauge and pseudofermion fields. A fictitious mo- 



mentum is introduced for each link variable with a Gaus- 
sian distribution, and a Hamiltonian is formed involving 
these momenta and the original action as a potential en- 
ergy. A new field configuration is proposed by approxi- 
mately solving Hamilton's equations for some length of 
fictitious time, known as an RHMC trajectory. In each 
ensemble, successive configurations are separated by 20 
RHMC trajectories to minimize autocorrelations. An 
improved anisotropic clover fermion action and an im- 
proved gauge field action are used [111 . In these ensem- 
bles, /3 = 1.5 and the s quark mass parameter is set to 
TOs — —0.0743 in order to reproduce a specific combi- 
nation of hadron masses [28j. In the (24'^|390) ensemble, 
the light quark mass parameters are set to m,j = = 
—0.0840 so that the pion mass is around 390 MeV if one 
sets the scale using the Vl baryon mass. In the (32'^|240) 
and (24'^|240) ensembles, m,j = md — —0.0860 are used, 
resulting in a pion mass around 240 MeV. The spatial 
grid size is ^ 0.12 fm, whereas the temporal spacing 
is at ^ 0.035 fm. 

A description of step (b) above, computation of the 
quark sinks, is given in Ref. Q. We employ the Lapla- 
cian Heaviside (LapH) quark-field smearing scheme de- 
fined using the three-dimensional gauge-covariant Lapla- 
cian expressed in terms of a stout-smeared gauge field. 
The spatial links are smeared using the stout-link proce- 
dure described in Ref. [13] with = 10 iterations and 
staple weight ^ = 0.10. For the cutoff in the LapH smear- 
ing, we use fj^ = 0.33, which translates into the number 
of LapH eigenvectors retained being — 112 for the 
24^ lattices and 7V„ = 264 for the 32^ lattice. We use Z4 
noise in all of our stochastic estimates of quark propa- 
gation. Our variance reduction procedure is similar to 
that described in Ref. [l^l- Our noise dilution projectors 
are products of time dilution, spin dilution, and LapH 
eigenvector dilution projectors. We use a triplet (T, S, 
L) to specify a given dilution scheme, where "T" denotes 
time, "S" denotes spin, and "L" denotes LapH eigenvec- 
tor dilution. The schemes are denoted by 1 for no dilu- 
tion, F for full dilution, and BK and IK for block-X and 
interlace- A', respectively (see Ref. |6|]). For all forward- 
time quark lines connecting source time to to the later 
sink time tp, we use the dilution scheme (TF, SF, Llj), 
where j — 8 for mesons and j = A for baryons. We 
chose J = 4 for baryons in order to dramatically reduce 
the disk space needed to store the baryon sources/sinks 
by a factor of eight. We found that the statistical er- 
rors in the baryon effective masses only increased by a 
factor of two in changing from LI8 to LI4. Given the 
eightfold reduction in both storage and computing time, 
we deemed this an acceptable loss of accuracy. For all 
same-sink-time tp-io-tp quark lines, the dilution scheme 
(TI16, SF, Llj) is used, where j — 8 for mesons and j = 4 
for baryons. The TI16 interlacing in time enables us to 
evaluate quark lines that originate on any time slice, al- 
lowing us to evaluate all diagrams needed to obtain the 
temporal correlations involving single-hadron and multi- 
hadron operators. Four widely-separated source times io 
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TABLE XX: The numbers A'op of single-hadron operators 
whose sinks and sources have been computed. The opera- 
tor numbers include operators having zero, on-axis, planar- 
diagonal, and cubic-diagonal momenta, as well as momenta 
for the pions in various special directions, such as (0, 1, 2) and 
(1,1,2). A'mom indicates the number of momenta in the op- 
erator sets. Here, tt refers to any isovector quark-antiquark 
meson operator with flavor content such as du (such as a, b, 
TT, and p mesons), rj refers to any uu + dd isoscalar meson op- 
erator, <f) refers to any ss isoscalar meson operator, K refers 
to any quark-antiquark operator having flavor content su or 
sd such that its strangeness is S — 1, and K is any quark- 
antiquark operator having flavor content us or ds such that 
its strangeness is S — —1. 





(24^1390) 


(24^1240) 


(32^1240) 




A^op 


mom 


A^op 


mom 


A^op 
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TT 
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2028 
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51 
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51 
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51 
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51 
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77 


K 
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51 
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51 
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65 


'K 
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51 


1517 


51 


1949 


65 


N/A 


1472 


33 






1616 


59 


A/E 


2274 


33 






2054 


51 




1320 


33 






700 


51 


Q. 


728 


33 






680 


45 



are used on each 2A^ gauge configuration, whereas eight 
to values are used on the 32'^ lattice. 

Details about how the single-hadron sources and sinks 
are computed in step (c) above can be found in Ref . [6] . A 
comprehensive survey of the spectrum of excited states in 
QCD requires obtaining temporal correlations of a large 
number of different operators. For example, KK oper- 
ators must mix with tttt operators, and so on. Hence, 
a large variety of single-hadron operators must be avail- 
able in the study of each sector, and software to carry 
out the plethora of different Wick contractions is needed. 
In order to study all stationary states of QCD involv- 
ing the u, d, s quarks, we need suitable sources and sinks 
for all isovector mesons (tt, p, a, 6), all isoscalar mesons 
{rj, f,h,u!,(l)), and all kaons. To study KK states, we 
also need separate antikaon sources and sinks. For the 
isoscalar mesons, separate uu + dd and ss operators 
are evaluated, and these will be allowed to mix in the 
correlation matrices. We also need all single-baryon 
A'', A, S, A, S, sources and sinks. Each hadron source 
and sink involves a summation over color indices and the 
spatial sites on each time-slice of the lattice. Different 
spin components and displacement directions must be 
combined to form the hadron operators which transform 
irreducibly under the symmetry operations of the spatial 
lattice. 

Using the lowest-lying energies in each symmetry sec- 
tor as determined from low-statistics runs on small 16'^ 
lattices, we systematically identified all momenta in the 



different little group irreps that would be needed to cap- 
ture the energy spectrum up to 0.5a^^, where at is the 
temporal lattice spacing. This energy value ensures that 
a sufficient number of hadrons accessible to experiments 
can be studied without overwhelming the computational 
resources available to us. The numbers A'op of single- 
hadron operators whose sinks and sources have been 
computed are listed in Table IXXl The operator numbers 
include operators having zero, on-axis, planar-diagonal, 
and cubic-diagonal momenta, as well as momenta for 
pions in various special directions, such as (0,1,2) and 
(1,1,2). iVmom indicates the number of momenta in 
the operator sets. These operators were chosen from 
much larger sets. Small-volume, low-statistics runs were 
done to select operators having lower statistical noise and 
smaller contamination from higher-lying eigenstates. In 
this table, n refers to any isovector quark-antiquark me- 
son operator with fiavor content such as du (such as a, 6, 
TT, and p mesons), 77 refers to any uu + dd isoscalar meson 
operator, refers to any ss isoscalar meson operator, K 
refers to any quark-antiquark operator having flavor con- 
tent su or sd such that its strangeness is = 1, and K 
is any quark-antiquark operator having flavor content Its 
or ds such that its strangeness is 5 = — 1. The N and A 
were computed simultaneously, and so were the A and E, 
since these baryons share many of the same three-quark 
components. 

Our final correlation matrices will not use such large 
numbers of operators. We will combine the single-hadron 
operators into two- and three-particle operators, pro- 
jected into the irreps of the Oh point group. A last 
round of operator selections will then occur. Many of 
the multi-hadron operators will not produce signals of 
suitable quality or be linearly dependent with other oper- 
ators. We will have to identify a final set of multi-hadron 
operators that best allows us to extract the low-lying en- 
ergy spectrum in each symmetry sector. 

For baryons containing identical quark flavors, the 
sources/sinks defined by Eq. (23) in Ref. [6] can be evalu- 
ated by assigning the quark spin and displacement indices 
to the identical-fiavor quark lines using either a single 
canonical one-to-one mapping or by averaging over all 
such one-to-one mappings. Averaging over the quark- 
line mappings is tantamount to averaging over different 
noises, so a reduction in the correlator variances can oc- 
cur. However, operators of this form have more three- 
quark terms in them, and the resulting increase in the 
number of "elemental" three-quark operators that must 
be evaluated raises the computational costs. We did some 
small lattice studies to compare the benefits/cost of the 
two different ways of assigning spin and displacement 
indices to the identical-fiavor quark lines. We decided 
that the modest increase in computational cost was war- 
ranted, given the variance reduction that resulted. 

Since our hadrons are made out of spatially-displaced 
quark fields, we carried out a systematic study of the 
role of the displacement length in the meson and baryon 
operators. Excited-state contamination in the correla- 
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tors of various meson and baryon operators was investi- 
gated for a range of displacement lengths. For the value 
of fls ^ 0.12 fm in the ensembles used, a displacement 
length of Sfls produced operators having better overlaps 
with the low-lying mesons of interest, whereas displace- 
ment lengths 2as and 3as worked well for baryons. These 
displacements can be rather costly for baryons, so we de- 
cided to use 2as for all baryon operators. 

An outline of how the temporal correlations are ob- 
tained from the hadron sources and sinks, step (d) above, 
is given in Ref. [Q]. Evaluating the thousands of single- 
hadron operators having different momenta for different 
noises, dilution schemes, and types of quark line ends 
(source or sink, normal mode or 75 Hermitian conjugate 
mode using Eq. (17) in Ref. [6[) is done out of necessity in 
multiple separate runs on XSEDE computing resources, 
resulting in the sources and sinks of each operator being 
distributed across a multitude of different files. The first 
task in step (d) involves data reorganization to collect 
together the different noises and quark line ends for each 
needed single hadron into a single file. The multi-hadron 
sources and sinks are also constructed in this first task of 
step (d) . The coefficients to combine the different single- 
hadron lines into multi-hadron operators are read from 
file and used to build and store the needed sources and 
sinks for each multi-hadron operator in a single file. A 
special data format is used that dramatically reduces the 
storage space needed for the multi-hadron sources and 
sinks. 

The second task in step (d) is computing the tempo- 
ral correlations using expressions such as Eqs. (24) and 
(33) in Ref. Q. Given the large number of correlators 
that we will need to evaluate, it is important to auto- 
mate the Wick contraction process as much as possible 
and to develop efficient software to handle all isospin sit- 
uations. Maple code was developed to create the C-|— I- 
subroutines for use in the numerical evaluations. The 
first step in the Maple calculations is forming appropri- 
ate flavor combinations of symbolic Grassmann variables 
representing single-hadron and multi-hadron operators of 
definite total isospin /, strangeness S*, and isospin pro- 
jection /a, raising and lowering the isospin projections 
of the individual constituent hadrons as needed. Next, 
the Maple code carries out the actual Wick contractions 
in terms of the Grassmann symbols. Finally, our Maple 
program outputs the C-f--|- subroutines for use in numer- 
ically evaluating the correlators in terms of the stored 
hadron sources and sinks in the different flavor sectors. 
For each flavor sector, we obtain a subroutine that eval- 
uates an expression similar to those shown in Eqs. (24) 
and (33) in Ref. @. Development and testing of the re- 
sulting C-l— I- software suite, known as "lastJaph", has 
now been completed. In evaluating the individual quark- 
line diagrams, we apply 75 Hermiticity (see Eq. (16) in 
Ref. 0) in cases where a V'(io) at the source connects 
with a ipitp) at the sink. For all same-time quark lines, 
we average over estimates obtained using both Eqs. (15) 
and (16) in Ref. 6] to increase statistics. 



IV. TESTS OF THE SINGLE-HADRON 
OPERATORS 

In order to test the effectiveness of the single-hadron 
operators that we have designed, we examined the effec- 
tive masses associated with the correlators of a variety 
of meson and baryon operators having various momenta. 
Samples of these effective masses are shown in Figs. [2] 
and [3] In these figures, we use the following definition of 
the effective mass: 

where time separation t in the correlator C (t) is measured 
in term of the temporal lattice spacing at, and usually 
dt ~ 3 is used. 

Fig. [2]shows the effective masses, mes{t) using dt — 3, 
associated with several isovector and isoscalar meson op- 
erators of various momenta p = {2tt/L) piat on the 
(24'^|390) ensemble, where L — 24as is the spatial ex- 
tent of the lattice. SS refers to a single-site meson oper- 
ator, and LSD denotes an operator in which the quark 
is displaced from the antiquark in a longitudinal direc- 
tion along the direction of the momentum, as explained 
in Sec. |lll In a TSD operator, the quark is displaced in 
a direction transverse to that of the momentum. Sev- 
eral different operators of each spatial type can be con- 
structed, and we label these different operators using an 
integer that varies from zero to the number of such op- 
erators less one. This integer is placed at the end of 
the spatial type label, such as SSO or TSD4. Results for 
on-axis momentum piat = (0,0,1), planar-diagonal mo- 
mentum Plat = (0,1,1), and cubic-diagonal momentum 
Plat = (1,1,1) are shown. The names of the lowest- lying 
known particles which appear in these finite- volume sym- 
metry channels are indicated in square brackets in each 
plot. 

Energies of the pion, nucleon, and A, obtained using 
the (32^ 1 240) ensemble, are studied in Fig. |31 In the left 
plot of this figure, effective masses associated with single- 
site pion operators having different momenta are shown 
against temporal separation t using dt = 3. Each effec- 
tive mass is labelled by its piat. The dispersion relation 
for the pion is shown in the center plot, which displays 
a^E'^ against pf^^, where E is the pion energy. Effective 
masses associated with two single-site baryon operators 
having zero momentum are shown in the right plot of 
Fig. [3] One effective mass corresponds to the nucleon 
(Gig channel) and the other to a A baryon {Hg chan- 
nel). 

Both Figs. [2] and [3] show that the stochastic LapH 
method works well for moving hadrons on large lattices. 
These plots also confirm that our choices of the gauge- 
field and quark-field smearing parameters and our noise 
dilution schemes for reducing variances are suitable for 
treating hadron states having nonzero definite momenta 
on the three Monte Carlo ensembles we plan to use. 
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FIG. 2: Effective masses, meff(t) using dt — 3, associated with several isovector and isoscalar meson operators of various 
momenta p — {2-jt/L) piat on the (24'^|390) ensemble, where L — 24as is the spatial extent of the lattice. SS refers to a 
single-site meson operator, and LSD denotes an operator in which the quark is displaced from the antiquark in a longitudinal 
direction along the direction of the momentum, as explained in Sec.|ll] In a TSD operator, the quark is displaced in a direction 
transverse to that of the momentum. The numbers following the letters, such as in SSO, are simply identifying integer indices. 
Particle names in the square brackets have been included to lend context to the little group irrep labels. 



TWO-HADRON OPERATORS OF DEFINITE 
MOMENTUM 



We construct our two-hadron operators using the same 
procedure described in Sec.|TT]to build the single-hadron 
operators, except that the basic building blocks are 
now the single-hadron operators instead of the elemental 
three-quark and quark-antiquark operators. 

Each single-hadron operator is labelled by total isospin 



/, the projection of the total isospin I^, strangeness S, 
three-momentum p, the little group irrep A, the row of 
the irrep A, and i, which denotes all other identifying 
information, such as the displacement type and index. 
Hence, basis operators for the two-hadron operators can 
be written 



(28) 



where B denotes either a baryon (as a Grassmann opera- 
tor) or a meson, and a and b denote the separate hadrons. 
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FIG. 3: Pion and baryon results obtained using the (32"^j240) ensemble. (Left) Effective masses, ma«{t) using dt = 3, against 
temporal separation t, associated with single-site pion operators having different momenta. Each effective mass is labelled by its 
Plat. (Center) Dispersion relation for the pion showing a^E"^ against pf^ti where E is the pion energy. (Right) Effective masses, 
mefi{i) using dt = 3, against temporal separation t, associated with two single-site baryon operators having zero momentum. 
One effective mass corresponds to the nucleon (Gig channel) and the other to a A baryon {Hg channel). 
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FIG. 4: Effective masses, mcfi(t) using dt — 3, associated 
with a two-meson operator in the Tiu irrep, having total 
isospin I = ^ and zero total momentum, constructed from 
single-site kaon and pion operators having equal and opposite 
on-axis momenta of minimal nonzero magnitude. Results on 
the (24^1390) and (32^ 1 240) ensembles are shown. 



Although the above operators form a perfectly acceptable 
basis of operators, they do not transform irreducibly un- 
der isospin rotations, nor under the 0\ symmetry trans- 
formations. 

To construct such irreducible operators, we first need 
to know how the basis operators transform. Eqs. (P^. 
((23)) . ((25)) . and (|26p summarize the important transfor- 
mation properties of our singic-hadron operators. Start- 
ing with the basis operators above, we first identify sub- 
sets of these operators that transform among themselves 



under all transformations. The total momentum 
p = Pa+Pb^s the first quantity we consider. We identify 
the little group of transformations that leave p invariant, 
then our goal is to construct operators that transform ac- 
cording to the A irrep of that little group. Under group 
element i? of O^, the single particle operator B^'^J^^Xaia 
transforms into one that has possibly a different momen- 
tum Rpa and is a linear combination of the different rows 
of the Aa irrep. The following quantities do not change 
under the transformations: Aa, Ab, ia, ib- We fix the 
above quantities, then apply the same group-theoretical 
projections as for the single-hadrons to construct the lin- 
ear combinations that transform irreducibly. This can 
be easily done since we know exactly how each of the ba- 
sis operators transforms under any 0\ transformation. 
Next, we form flavor combinations that transform ir- 
reducibly under isospin rotations. Lastly, we apply G- 
parity projections, whenever suitable. 

The above procedure is applied for two-particle opera- 
tors having total momenta in the directions of the refer- 
ence momenta Pref in Table |TT1 An operator having total 
momentum in any other direction is obtained by applying 
the appropriate reference rotation R^of in Table mi to the 
appropriate operator having momentum in a reference 
direction. 

In order to test the effectiveness of the two-hadron op- 
erators that we have designed, we examined the effective 
masses associated with the correlators of a variety of two- 
hadron operators. We also evaluated several correlation 
matrices mixing single and two-hadron operators. Sam- 
ples of our results are shown in Figs. 21 [5l and El 

In Fig.m effective masses using dt — 3 associated with 
a two-meson operator in the Tiu irrep are shown. The 
two-meson operator has total isospin I = ^ and zero 
total momentum and is constructed from single-site kaon 
and pion operators having equal and opposite on-axis 
momenta of minimal nonzero magnitude. Results on the 
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FIG. 5: Mixing of a quark-antiquarlc single-site p operator and a two-pion operator. Both operators transform according to 
the Tjl^ irrep and create states of zero total momentum. In the two-pion operator, each pion is a single-site operator travelling 
with minimal nonzero on-ajcis momenta. (Left) The 2x2 correlation matrix Cij{t) ( Cii{ro)Cjj{ro) )~^/^ with tq = 3 in terms 
of temporal lattice spacing at. (Center) Effective masses, mcft(t) using dt — 3, associated with the diagonal elements of the 
original correlator matrix. (Right) Effective masses associated with the optimized correlator (see text). The meanings of the 
horizontal solid and dashed lines are described in the text. These results were obtained on the (24^1390) ensemble. 



(243|390) and (323|240) ensembles are shown. 

The mixing of a single-site p quark-antiquark opera- 
tor and a two-pion operator is examined in Fig. [5j Both 
operators transform according to the T^^ irrep and cre- 
ate states of zero total momentum. In the two-pion 
operator, each pion is a single-site operator travelling 
with minimal nonzero on-axis momenta. In the left 
plot, the 2x2 real and symmetric correlation matrix 
Cij{t) ( Cu{To)Cjj{To) )~^/^, with To = 3, is shown. Ir- 
relevant normalization factors are removed by dividing 
the original correlation matrix by ( Cii{To)Cjj{To) Y^'^ 
for some early tq. In each matrix element label, the sym- 
bol to the right of the vertical bar indicates the source 
operator and the symbol to the left of the vertical bar 
denotes the sink operator. These results were obtained 
on the (24'^ 1 390) ensemble. This figure provides further 
evidence that the stochastic LapH method is well suited 
to studying systems involving the mixing of single and 
multi-hadron operators, and that our smearing parame- 
ters and noise dilution schemes were chosen judiciously. 
One sees that there is a small but noticeably nonzero 
mixing between these operators. 

In the center plot of Fig. [31 the effective masses associ- 
ated with the diagonal elements of the original correlator 
matrix are displayed. A so-called "optimized" correla- 
tion matrix C°P*(t) can be defined by first determining 
the unitary matrix U in the eigenvalue equation below: 

C(to)-i/2 C'iTD) C(to)-i/2 = u C°P'(rc) C/t, (29) 

where C°p*(t£)) is a diagonal matrix and C(t) is the orig- 
inal Hermitian correlation matrix, then writing 

C°P*(t) ^ U^C'iTo)-^/^ C{t) C(to)-i/2 U. (30) 

Here, tq = 8 and td = 15 were used to ensure that the 
optimized correlation matrix remains diagonal for all t 



exceeding td. In the right plot in Fig. [SJ the effective 
masses associated with the diagonal elements of C°P*(i) 
are shown. 

The main purpose of Fig. [S] is to demonstrate the ef- 
fectiveness of the stochastic LapH method in providing 
estimates of correlation matrices involving both single 
and two-hadron operators that are accurate enough to 
allow reliable diagonalizations and extractions of excited- 
state energies. However, the center and rightmost plots of 
Fig. [5] also illustrate a major pitfall that must be avoided 
when trying to extract excited-state energies in lattice 
QCD calculations. The effective masses of the optimized 
correlators tend to the two lowest-lying stationary-state 
energies in this symmetry channel, as they should. The 
solid horizontal lines, shown in both plots, indicate the 
location of the first excited-state energy as determined 
by a single-exponential fit, with wrap-around, to the op- 
timized correlator. An excellent fit quality is obtained, 
and the uncertainty in the energy from the fit is indi- 
cated by the two parallel lines. In the center plot, one 
sees that the effective mass associated with the original 
p operator tends to the lowest-lying stationary-state en- 
ergy in this channel, consistent with expectations. The 
pitfall is revealed in the effective mass associated with 
the two-pion operator in the center plot. This effec- 
tive mass deceptively appears to level off at a plateau, 
shown by the dashed lines in both the center and right- 
most plots. The dashed lines show the result of fitting 
a single-exponential fit, with wrap-around, to this cor- 
relator, and a good quality of fit is obtained. However, 
this operator has a small but nonzero coupling to the 
lowest-lying state, so the dashed line must be viewed as 
a "false plateau" . The solid line indicates the actual en- 
ergy of the first-excited stationary state in this channel, 
as determined from the diagonalized correlation matrix, 
which lies above this false plateau in the effective mass. 
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FIG. 6: The 5x5 matrix Cij{t) ( C«(ro)Cjj(ro) )~^^^, with to = 4, of the temporal correlations of five isoscalar operators in 
the scalar A^^ sector for the (24"^|390) ensemble. G is the Ga glueball operator described in Sec. lVIl and a refers to a single-site 
quark-antiquark operator. Three of the operators are 7 = two-meson operators constructed out of single-site single-meson 
operators having equal and opposite momenta. In the r]{0)ri{0) and 7r(0)7r(0) operators, each meson has zero momentum, 
whereas each pion in the 7r(l)7r(— 1) has minimal nonzero on-axis momentum. Note that large vacuum expectation values have 
been subtracted to obtain each of these correlators. 



This shows that it is unwise to ignore small couplings 
to lower-lying states. The effects of these small couplings 
are usually not different enough on different lattice vol- 
umes to allow their recognition. This serves as a strong 
warning that to extract a resonance lying above a multi- 
particle threshold, complete removal of overlaps onto the 
lower-lying multi-hadron states is mandatory, and this 
usually requires the introduction of good multi-hadron 
operators into the correlation matrix analysis. This issue 
will be investigated in more detail in future work. 

A last example of the effectiveness of the stochastic 
LapH method in dealing with mixings between single and 
two-meson operators is given in Fig. [51 which shows a 
5 X 5 correlation matrix Cij{t) ( Cii(ro)Cjj (tq) )~^^^, with 
To = 4, of the temporal correlations of five isoscalar oper- 
ators in the scalar A^^ sector for the (24'^ 1 390) ensemble. 
G is the Ga glueball operator described in Sec. lVIl and a 
refers to a single-site quark-antiquark operator. Three of 



the operators are 7 = two-meson operators constructed 
out of single-site single-meson operators having equal and 
opposite momenta. In the 77(0)77(0) and 7r(0)7r(0) op- 
erators, each meson has zero momentum, whereas each 
pion in the 7r(l)7r(— 1) has minimal nonzero on-axis mo- 
mentum. Note that large vacuum expectation values 
have been subtracted to obtain each of these correlators. 
This is a notoriously difficult channel to study, but the 
stochastic LapH method appears to produce results of ad- 
equate precision, even for the rjrj — rji] correlator, which 
includes many diagrams involving internal quark loops 
and other same-time quark lines. The a operator does 
not appear to mix much with the other four operators. 
We certainly plan to investigate this channel in much 
more detail in the future. 

Table IXXll summarizes the different fiavor types of sin- 
gle and two-hadron operators that we plan to include in 
our first survey of the spectrum of stationary-state en- 
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TABLE XXI: The following bosonic and fermionic flavor sec- 
tors involving only the u, d, s quarks will be studied. J de- 
notes total spin, / denotes total isospin, and S is the total 
strangeness. The particle contents to be studied in each fla- 
vor sector are listed. G denotes a glueball operator. Ifere, vr 
refers to any isovector quark-antiquark meson operator with 
flavor content such as du (such as a, b, n, and p mesons), 
77 refers to any uu + dd isoscalar meson operator, </> refers 
to any ss isoscalar meson operator, K refers to any quark- 
antiquark operator having flavor content su or sd such that 
its strangeness is S = 1, and K is any quark-antiquark oper- 
ator having flavor content us or ds such that its strangeness 
is 5* = -1. 
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n to 


study all bosonic and fermionic fla- 



vor sectors involving the u, d, s quarks that involve up to 
two meson and meson-baryon pairs. The particle types 
to be studied in each flavor sector are listed in this ta- 
ble. G denotes a glueball operator. In this table, tt 
refers to any isovector quark-antiquark meson operator 
with flavor content such as du (such as a, b, tt, and p 
mesons), 77 refers to any uu + dd isoscalar meson opera- 
tor, (f) refers to any ss isoscalar meson operator, K refers 
to any quark-antiquark operator having flavor content su 
or 'sd such that its strangeness is S" = 1, and K is any 
quark-antiquark operator having flavor content us or ds 
such that its strangeness is = — 1. Our current plans 
do not include states containing two or more baryons or 
three or more mesons, although the computational tech- 
nology can easily accommodate such states. 

VI. A NEW GLUEBALL OPERATOR 

Determining stationary-state energies in the interest- 
ing scalar isoscalar sector will ultimately involve includ- 
ing a scalar glueball operator. Glueballs are hypotheti- 



cal particles comprised predominantly of gluons, having 
no valence quarks. Scalar glueball operators are usually 
constructed using a sum of gauge-invariant loops of the 
smeared spatial link variables on a single time slice, which 
is invariant under translations, rotations, and charge con- 
jugation. However, any purely gluonic quantity with 
similar symmetry properties could presumably be used. 
LapH quark-field smearing involves the covariant spatial 
Laplacian A. The eigenvalues of the Laplacian are in- 
variant under rotations and gauge transformations, and 
so, are appropriate for a scalar glueball operator. The 
lowest-lying eigenvalue was studied, as well as other func- 
tions of the eigenvalues. We found that essentially any 
combination of the low-lying eigenvalues worked equally 
well for studying the scalar glueball. Two operators in 
particular that we studied are 

GaW = -Tr[e(a2-f A) A], ^ (31) 

Gw{t) = -Tr[e(a2 + A) A exp(-M^A2)]. (32) 

The first operator Ga in Eq. (I5l]) . which we call the 
TrLapH operator, is perhaps the simplest operator that 
one can construct using the eigenvalues of the covariant 
Laplacian. In the so-called weighted TrLapH operator 
Gw in Eq. (P^ . we used = 64 in order that only a 
handful of the lowest-lying eigenvalues contribute. 

The effective masses obtained using Ga and Gw on 
the 24^ X 128 ensemble are shown in Fig. [71 These effec- 
tive masses are compared to that obtained using a stan- 
dard glueball operator which is a sum of 3 x 3 loops of 
the smeared gauge link variables that is rotationally and 
translationally invariant. One observes very little dif- 
ference between these effective masses, suggesting these 
operators are comparable in usefulness for studying the 
scalar glueball. Similar conclusions were reached using 
the 32^^ X 256 lattice. Each effective mass eventually 
tends toward the energy of two pions at rest, demon- 
strating non-negligible coupling of these operators to tttt 
states. Thus, we plan to use the simplest operator G^it) 
in future studies involving the scalar glueball. 

VII. CONCLUSION 

Multi-hadron operators are crucial for reliably ex- 
tracting the masses of excited states lying above multi- 
hadron thresholds in lattice QCD Monte Carlo calcula- 
tions. Multi-hadron operators with significant coupling 
to the low-lying states of interest can be obtained by 
combining single-hadron operators of various momenta. 
The construction and testing of single-hadron operators 
of definite momentum, and their combinations into two- 
hadron operators was the main subject of this work. 

The approach of Ref. [l| was extended to meson oper- 
ators of zero momentum, and to both meson and baryon 
operators having definite nonzero momentum. Our op- 
erator design utilizes group-theoretical projections. The 
point and space groups we use are well known, and the 
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FIG. 7: Comparison of the effective masses, mcff(f) using dt = 3, associated with three different scalar glueball operators on 
the (24''|390) ensemble, (a) The leftmost plot shows the effective mass using an operator defined by a sum of 3 x 3 loops of 
the smeared gauge link variables that is rotationally and translationally invariant, (b) The middle plot uses the new TrLapIf 
glueball operator Ga, defined in Eq. (|3ip . (c) The rightmost plot shows the results using the weighted TrLapH glueball operator 
Gw, defined in Eq. (|32|l . One observes very little difi^erence between these plots, suggesting that these operators are comparable 
in usefulness for studying the scalar glueball. Each effective mass eventually tends toward the energy of two pions at rest, 
demonstrating non-negligible coupling of these operators to tttv states. 



properties of their irreducible representations are widely 
available in the literature. However, we collected to- 
gether and presented in this paper some of the specific 
group theory details needed for our operator construction 
for the convenience of the reader and to explicitly state 
our conventions and the notations we use. 

Tests of our single-hadron operators using a stochastic 
method of treating the low-lying modes of quark prop- 
agation which exploits Laplacian Heaviside quark-field 
smearing were presented. These tests were carried out 
on 24^ X 128 and 32'^ x 256 anisotropic lattices with pion 
masses m-^ w 390 and 240 MeV. These tests revealed 
a pitfall that must be avoided when trying to extract 
excited-state energies in lattice QCD, namely, that the 
small couplings of interpolating operators to lower-lying 
states must not be ignored. This issue will be further 
studied in future work. A new glueball operator was also 
introduced and tested. We demonstrated that comput- 
ing the mixing of this glueball operator with a quark- 
antiquark operator, tttt, and rjri operators is feasible with 
the stochastic LapH method. 

The stochastic LapH method provides reliable esti- 
mates of all temporal correlations that will be needed 
for a comprehensive survey of the low-lying spectrum of 
QCD stationary states in finite volume. The method 
works well even for those correlators that are particu- 
larly difficult to compute, such as rjri — 7777 in the scalar 
channel, which involves the subtraction of a large vac- 
uum expectation value. The effectiveness of the method 
can be traced to two of its key features: the use of noise 
dilution projectors that interlace in time, and the use of 
Zjv noise in the subspace defined by the Laplacian Heavi- 
side quark-field smearing. Introducing noise in the LapH 
subspace results in greatly reduced variances in temporal 
correlations compared to methods that introduce noise 
on the entire lattice. Although the number of Laplacian 



eigenvectors needed to span the LapH subspace rises dra- 
matically with the spatial volume, the number of inver- 
sions of the Dirac matrix needed for a target accuracy 
is remarkably insensitive to the lattice volume, once a 
sufficient number of dilution projectors is introduced. 

In addition to increased efficiency, the stochastic LapH 
method has other advantages. The method leads to com- 
plete factorization of hadron sources and sinks in tempo- 
ral correlations, which greatly simplifies the logistics of 
evaluating correlation matrices involving large numbers 
of operators. Implementing the Wick contractions of the 
quark lines is also straightforward. Contributions from 
different Wick orderings within a class of quark-line di- 
agrams differ only by permutations of the noises at the 
source. 

In the future, we plan to carry out a comprehensive 
survey of the excitation spectrum of the stationary states 
of QCD involving mesons and baryons containing u, d, s 
quarks. Various scattering phase shifts and decay con- 
stants will also be investigated. The needed single-meson 
and single-baryon sources and sinks for a large number 
of different momenta have been computed and stored for 
three Monte Carlo ensembles. The development and test- 
ing of the software to combine these sources and sinks via 
Wick contractions of the quark fields into temporal cor- 
relators has been completed, and the final stages of our 
operator selections are now in progress. Results for the 
spectrum, involving both single- and two-hadron opera- 
tors and using the technology described in this work, will 
appear in future publications. 
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